Correlation between event multiplicity and elliptic flow parameter 
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The elliptic flow parameter (component), V2, in the Fourier expansion of event-by-event charged 
particle multiplicity azimuthal distribution in the momentum space is studied by taking into account 
the multiplicity fluctuations. The correlations between charge multiplicity and impact parameter 
as well as between elliptic flow parameter and impact parameter (multiplicity) are investigated by 
both the parton and hadron cascade model PACIAE and a multiple phase transport model AMPT. 
A discussion is given for the event-wise average and particle-wise average in the definition of elliptic 
flow parameter. 

PACS numbers: 25.75.Dw, 24.85.+p 
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I. INTRODUCTION 



The charge particles emitted from the fireball created 
in relativistic nucleus-nucleus collisions exhibit trans- 
verse collective flow. This is represented by the elliptic 
flow parameter (i^) and other harmonic parameters (v n , 
n=0, 1, 3, 4, ...). Those parameters, as Fourier expansion 
coefficients of produced particle azimuthal distribution in 
the momentum space, are highly sensitive to the spatial 
geometry (eccentricity) of the created fireball (nuclear 
overlap region or interaction region). 

The expected phase transition to Quark-Gluon-Plasma 
(QGP) should have a dramatic effect on those harmonic 
parameters. The consistency between experimental data 
of V2(pr) ( V2(y)) at mid-rapidity and the corresponding 
hydrodynamic predictions is regarded as an evidence of 
the production of partonic matter in the ultra-relativistic 
nucleus-nucleus collisions The elliptic flow of high 

Pt particles may be related to jet fragmentation and par- 
ton energy loss [3J] , which are usually not included in the 
hydrodynamic calculations. Such kind of hydrodynamic 
calculation Q overestimates V2(pt) in the pr > f.5 
GeV/c region This is regarded as an evidence of 
the strongly coupled QGP formation in the relativistic 
nucleus-nucleus collisions together with the discovery of 
jet quenching Q. So far a lot of experimental data have 
been published on the collective flow parameters @, B @ • 
Consequently microscopic transport model studies are 
also widely progressing [T(| El, E3, EH EH as well as 
the abundant hydrodynamic investigations. 

According to two well known pioneering works in this 
field " 
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the usual study starts from the triple dif- 
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[1+ 2« n cos[n(0-* r )]], (f) 



where TV is the particle multiplicity distribution, <fi stands 
for the azimuthal angle of particle, and \I/ r refers to 
the azimuthal angle of reaction plane in the momen- 
tum space. Then the n-th flow harmonics is defined 
as v n = (cos[n((j) — 'J'r-)]), where () indicates an av- 
erage over all particles in all events of a sample. For 
the distribution of all particles in the sample, the coeffi- 
cients in Eq. Jl} can be calculated as v\ = (p x /pt) and 
v 2 = ((Px/pt) 2 — (Py/pr) 2 ), etc. This kind of average 
is widely accepted and will be indicated as particle-wise 
average hereafter. 

The problem of the particle-wise average in a sample 
with wider multiplicity range is that it does not take the 
influence of multiplicity (hence impact parameter) in a 
single event of the sample into account. In [l7| the elliptic 
flow and other harmonics have been re-derived starting 
from the invariant particle multiplicity distribution in the 
momentum space. It turned out that the harmonic pa- 
rameters v n (elliptic flow parameter V2) is an event- wise 
average of cos(n<p) 



v e n = (cos(n<t>)) ev (n = 1,2, ...), 



(2) 



where cos(n4>) denotes the average of cos(n(j)) over parti- 
cles in a single event, (...) &v means an average over events 
in a sample, and the superscript "e" stands for the event- 
wise average. It is an average of cos(n<t>) first over par- 
ticles in a single event then over the events in a sample. 
Here it has to mention that in theory if the beam direc- 
tion and impact parameter vectors are fixed at the p z and 
p x axes, respectively, then the reaction plane is just the 
Px ~Pz plane (l5j . Therefore the reaction plane azimuthal 
angle in Eq. (JTJ) introduced for the extraction of el- 
liptic flow in experiments [l6| is zero. Meanwhile, the 



2 



particle-wise average of v n is also derived in [lj 



nP = 



cos(n<p)N ev ) I (N ev ) e 

I ev 



(3) 



That is obviously different from event- wise average. Only 
if the cos(n(f>) is independent of event multiplicity N ev (i. 
e. if the multiplicity plays no role in the average) the 
reduces to v„. In fact, the cos(n(f>) and N ev correlate 
(even negatively correlate) with each other. This is be- 
cause larger event multiplicity arises from more central 
collisions (larger overlap region between colliding nuclei) 
and the larger overlap region, in turn, results in less az- 
imuthal asymmetry. The particle-wise average does not 
take the influence of event multiplicity into account, thus 
it is questionable from physics point of view. Of course, 
for a very narrow multiplicity bin studied, the particle- 
wise average is not very problematic. However for the 
wide multiplicity bin the correction is important. Here 
we do not take the influence of eccentricity fluctuation 
on i>2 into account, it may strengthen the importance of 
the correction. 

In experiment, the reaction plane is different event by 
event. In order to extract the elliptic flow parameter 
one has to invoke a complex reaction plane identifica- 
tion method the cumulant method [Hj|, and the 
Lee- Yang zeroes method [ijj]. In all of these methods a 
quantity has to be first constructed event by event. This 
quantity is just the event plane in [l6j . the cumulant ex- 
pansion of the weighted n-th transverse event- flow vector 
in (l8| . and a generating function in pj| . Then a corre- 
sponding average over measured events has to be taken. 
Therefore the experimental extraction of elliptic flow pa- 
rameter is an event-wise average. In the recent paper 
[20| it has been pointed out that "Elliptic flow develops 
on an event-by-event basis. However, the experimental 
determination of V2 demands averaging over events and 
thus over distributions of geometric shapes." This means 
again the experimental extraction of v n is an event-wise 
average. 

In the expansion process of a nucleus-nucleus colli- 
sion, the initial spatial asymmetry of the created fireball 
evolves into an azimuthal asymmetry (v n ) in the trans- 
verse momentum distribution of produced charged parti- 
cles This evolution process has to be described pre- 
cisely by detailed dynamical models. Thus, in this work 
a parton and hadron cascade model PACIAE [2l| and a 
multiple phase transport model AMPT (with string melt- 
ing) [22j are both used in order to have cross checking. 
The observables of the impact parameter (6), charged 
particle multiplicity (N c h), eccentricity (e), and elliptic 
flow parameter (V2) as well as other harmonic parame- 
ters (v n ) are always used to describe the physics of the 
initial spatial asymmetry, final momentum asymmetry, 
and their fluctuations and correlations. As first step, in 
this paper we only study the correlation between N c h 
and b as well as the correlation between vi and b (N c h). 
In addition, the event-wise average versus particle-wise 
average in calculating (extracting) the elliptic flow pa- 



rameter raised in [l7| is also studied. The eccentricity 
fluctuation may have extra effect on the vi parameter 
and its fluctuation besides the impact parameter (mul- 
tiplicity) fluctuation, which will be studied in next step. 
Our studies are based on the multiplicity fluctuation at 
fixed impact parameter, and would apply also if the ini- 
tial eccentricity would not fluctuate. 



= 200 GeV 



N ch ( |T|I <l,0<p T <2GeV)) . 



AMPT 



b [fm] 



5 10 
b [fm] 



FIG. 1: Correlation between impact parameter b and charge 
multiplicity N c h calculated by the PACIAE model (left panel) 
and the AMPT model (right panel) for the indicated pseudo- 
rapidity and pt range. The error bars indicate the fluctua- 
tions of multiplicity at fixed impact parameter. Consequently 
a fixed event multiplicity may correspond to different impact 
parameters. 



II. MODELS 

The PACIAE model [2l| is a parton and hadron cas- 
cade model and consists of four stages: parton initializa- 
tion, parton evolution (rescattering), hadronization, and 
hadron evolution (rescattering). 

In the first stage, a nucleus-nucleus collision is decom- 
posed into nucleon-nucleon (NN) collisions according to 
the collision geometry: 

• Both colliding nuclei settle down at the distance 
between two centers of colliding nuclei along the 
x-axis is equal to impact parameter b. 

• Nucleons in the colliding nucleus are randomly dis- 
tributed around the center of nucleus according to 
a Wood-Saxon distribution (r) and the 4-7T uniform 
distribution {0 and 4>). 

• The beam momentum is given to nucleons in the 
colliding nucleus. The nucleon Fermi motion is ne- 
glected. 

• The NN collision is sampled randomly according 
to the spatial and momentum distributions, the 
straight line trajectories, and the NN interaction 
cross sections of the colliding pairs. 
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Then a NN collision is described by the PYTHIA model 
[H[. In the PYTHIA model 

• The NN collision is decomposed into parton-parton 
collisions. 

• A hard parton-parton collision is described by 
the lowest leading order perturbative QCD (LO- 
pQCD) parton-parton cross sections modified by 
the parton distribution function in a nuclcon. 

• The soft parton-parton collision is considered em- 
pirically. 

• Because the initial- and final-state QCD radia- 
tions are considered in parton-parton interactions 
the consequence of a NN collision is partonic mul- 
tijet configuration composed of di-quarks (anti- 
diquarks), quarks (anti-quarks), gluons, and a few 
hadronic remnants. 

• The string fragmentation is performed. Then one 
obtains a final hadronic state for a NN collision 
(hence for a nucleus- nucleus collision). 

However, in the PACIAE model the above string frag- 
mentation is switched-off temporarily. The di-quarks 
(anti-diquarks) are broken randomly into quarks (anti- 
quarks). So, the consequence of a NN collision (hence 
a nucleus-nucleus collision) is a configuration of quarks, 
anti-quarks, gluons, and hadronic remnants. In addition, 
there are spectator nuclcons for a nucleus- nucleus colli- 
sion. 

The partonic rescattering stage follows the initializa- 
tion one. In this stage the rescattering among partons 
are considered by the 2 — > 2 LO-pQCD differential cross 
sections [Hj]. The six elastic and three inelastic parton- 
parton scatterings are involved with their own differential 
cross sections [24j. We use the Monte Carlo method to 
simulate the parton rescattering until the parton-parton 
collisions cease. 

In the hadronization stage the partons after rescat- 
tering are hadronized by either the string fragmentation 
scheme (Lund string fragmentation and/or i ndep endent 
fragmentation) t 23] or the coalescence model [2l| . 

The last stage of hadronic rescattering is described by 
the usual two-body collision method [25|]. As such, the 
hadronic matter is evolved until the hadron-hadron col- 
lision pairs are exhausted. 

The two models of PACIAE and AMPT with string 
melting are physically similar. They have similar model 
assumptions but have otherwise very little commons in 
detail. The main differences between them are as follows: 

• AMPT is based on HIJING 0] instead of PYTHIA 
in PACIAE. 

• In the AMPT model the hadrons from HIJING are 
fragmented into partons, while the spectator nucle- 
ons are kept survival. This parton initialization for 
a nucleus-nucleus collision is quite different from 



the one in the PACIAE model, where parton ini- 
tialization is realized by broking the strings before 
the string fragmentation. Thus, even HIJING itself 
is heavily based on PYTHIA, the initialized parton 
state is quite different between the PACIAE and 
AMPT models. There are much more partons in 
the AMPT model than in PACIAE. 

• The Zhang parton cascade (ZPC) model 26] is em- 
ployed to describe parton rescattering. In ZPC only 
elastic scatterings are considered with gg — > gg 
cross section instead of all elastic interaction cross 
sections. 

• Partons after rescattering are hadronized by the 
coalescence model only. 

• The dynamics of the consequent hadronic matter is 
described by a relativistic transport (ART) model 
(27j . In ART the cross sections of hadron-hardon 
collisions and hadronic resonances are considered in 
detail. Therefore the hadronic rescattering is con- 
sidered more carefully in AMPT than in PACIAE. 

As mentioned in the beginning, the AMPT model has 
been used successfully to describe the elliptic flow param- 
eter [13, [3 in Au+Au collisions at RHIC energy. The 
AMPT results given in this paper are calculated by the 
same code [l2[ with parameters adjusted to the V2 data. 
However, the PACIAE results are calculated with the 
default parameters. Thus the two results are not match- 
ing each other very well. Nevertheless, that is irrelevant 
because we aim to explore the important physics men- 
tioned above rather than to reproduce the experimental 
data. To demonstrate that physics and to estimate its 
quantitative importance we used two independent mod- 
els, both having a realistic description for multiplicity 
fluctuations at the fixed impact parameter, at hand [28[ 
to cross check the correlation between N c h and b and the 
correlation between V2 and b (N c h). 



III. RESULTS 

The correlation between impact parameter b (in fm) 
and charged multiplicity N c h from the PACIAE and 
AMPT calculations for Au+Au collisions at y / S7v r /v=200 
GeV are shown in Fig. [1] One sees in the left panel of 
Fig. [T] that N c h is negatively correlated with b. A unit 
of "fm" change in impact parameter results in more than 
100 charged particle change in multiplicity. An about 
20% increase in multiplicity corresponds to about 2 fm 
decrease in impact parameter. Similar conclusions can 
be drawn from the AMPT results in right panel of Fig. 

m 
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FIG. 2: Impact parameter dependence of integrated V2 cal- 
culated by the PACIAE model (left panel) and the AMPT 
model (right panel). The error bars indicate the random fluc- 
tuations of vy at fixed impact parameter. 
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FIG. 3: Charged particle multiplicity dependence of inte- 
grated V2 calculated by the PACIAE model (left panel) and 
the AMPT model (right panel). The error bars indicate the 
random fluctuations of V2 at fixed charged particle multiplic- 
ity. 

In Fig. [2] we give the integrated vi (event-wise aver- 
age) as a function of impact parameter b calculated by the 
PACIAE and AMPT models for the Au+Au collisions at 
\/snn— 200 GeV. We see in left panel of Fig. [2]that there 
is a peak appearing at b nearly equal to the radius of col- 
liding nucleus. That is a result of competition between 
the charge multiplicity and the impact parameter. In 
central Au+Au collision the charge multiplicity reaches 
maximum at zero impact parameter, where nuclear over- 
lap region is nearly symmetric thus v-i approaches zero 
[B[ . In the middle peripheral collision, although the 
charge multiplicity is going down, the V2 is large because 
of the strong asymmetry of nuclear overlap region (large 
impact parameter). In the extra peripheral collision, the 
asymmetry of nuclear overlap region is very strong, while 
the V2 is going down because the multiplicity is too low 
to generate pressure gradient. When the impact param- 
eter is around the radius of colliding nucleus, the charge 
multiplicity is not so small, and the asymmetry of nuclear 
overlap region is considerably strong, so vi approaches its 
maximum. Similar situations can be seen in the AMPT 
results shown in the right panel of Fig. [51 These results 



are consistent with the PHOBOS and PHENIX reports 
about the peak structure in vi as a function of the num- 
ber of participant nucleons (N part ) in Au+Au collisions 
at v / saT/v=200 GeV [1, H§|- Similarly, we give the in- 
tegrated V2 as a function of N c h calculated by the PA- 
CIAE and AMPT models for the Au+Au collisions at 
v ^7=200 GeV in Fig. [3 

In the left panel of Fig. [T]we know that a nearly 20% in- 
crease in multiplicity corresponds to about 2 fm decrease 
in impact parameter. This change in impact parameter, 
in turn, results in about 0.010 change in v 2 (see the sen- 
sitive region of 2 < b < 7 in the left panel of Fig. \2§ . The 
similar, even stronger, conclusions can be drawn from 
the AMPT results in the right panels of Fig. [1] and [5] 
Thus, both the PACIAE and AMPT models, incorporat- 
ing random multiplicity fluctuation, verify that a charge 
multiplicity bin, which includes e.g. 20% of maximum 
charged multiplicity, spans an impact parameter range 
of Ab ~ 2 fm. That results in a change of Ai>2 ~ 0.010 in 
the PACIAE calculations. Considering the maximum vi 
in the PACIAE model is just only 0.028, so the change 
of Ai>2 caused by the different impact parameters in the 
sample is about 30-40% of the maximum! This is a very 
significant change, which should not be underestimated! 

In [l?]] the PACIAE model has been used to calcu- 
late the charged hadron v^irf) in 0-40% most central 
Au+Au collisions at GeV by the method of 

event-wise and particle-wise averages separately. The 
Vn(ri) is about 10% larger than vS(ri). 

is smaller than 



cos(n<p)N ev 



This means 

(Nev) e 



and 



demonstrates the negative correlation between cos(ncf>) 
and N ev . Here we use the AMPT model [22j to calculate 
i>f (77) and Ug^) in 0-40% most central Au+Au collisions 
at ^/i^^OO GeV again. The AMPT results of 1^(77) 
are nearly 20% larger than v^rj), as shown in Fig. [4] 
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FIG. 4: V2 as a function of rj calculated by the AMPT model 
for — 40% most central Au+Au at ^/snn= 200 GeV, using 
the two methods of averaging described in the text. 
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IV. CONCLUSIONS 

In summary, we have used the parton and hadron 
cascade model PACIAE and a multiple phase transport 
model AMPT with string melting to investigate the ellip- 
tic flow parameter u 2 by taking into account the event-by- 
event charge multiplicity fluctuations and to cross check 
the physics of initial spatial asymmetry, final momentum 
asymmetry, and their fluctuations and correlations. The 
correlations between the charge multiplicity and impact 
parameter as well as between the elliptic flow parame- 
ter and impact parameter (charge multiplicity) are cal- 
culated. 

It turned out that the charge multiplicity is negatively 
correlated with impact parameter. The v<i as a func- 
tion of impact parameter first increases with increasing 
6, reaches its maximum at b close to the radius of collid- 
ing nucleus, and turns to decrease with increasing b at 
last. This is because the is determined by two driving 
forces, multiplicity and impact parameter. 

The averaging procedure in the definition of elliptic 



flow parameter is reexamined by the AMPT calculations 
for vf (77) and v\(r\) in 0-40% most central Au+Au col- 
lisions at v / s^7=200 GeV. The AMPT results of uf (77) 
are about 20% larger than v^r/), which are larger than 
the corresponding PACIAE results [TtJ • This emphasizes 
again the necessity of theoretically calculating the ellip- 
tic flow parameter V2 by the event-wise average, in order 
to be consistent with the experimental extraction of V2, 
which is really event-wise average 



Acknowledgments 

Financial support from NSFC, China (grants 
10705012, 10475032, 10605040, 10635020, 10975062,) 
RCN, Norway (grants 101937 and 171247), the 
Chinese-Norwegian RCN-NFSC collaboration grant, 
no. 10710886, and Programme of Introducing Talents of 
Discipline to Universitiesnder under Grant No. B08033 
are gratefully acknowledged. 



[1] M. Gyulassy, |arXiv:nucl-th /04030321 and in "Structure 
and dynamics of elementary matter: Proceedings of 
the NATO Advanced Study Institute on Structure and 
Dynamics of Elementary Matter, Camyuva-Kemer (An- 
talya), Turkey, 22 September-2 October 2003, ed. by W. 
Greiner. 

[2] L.P. Csernai and D. Rohrich, Phys. Lett. B 458 454 
(1999). 

[3] M. Gyulassy and X.-N. Wang, Comput. Phys. Commun. 

83, 307 (1994). 
[4] P. F. Kolb, P. Huovinen, U. W. Heiz, and H. Heiselberg, 

Phys. Lett. B 500, 232 (2001). 
[5] B. B. Back, et al., Phobos collaboration, Phys. Rev. C 

72, 051901(R) (2005). 

[6] W. A. Zajc, larXiv:0802.3552l 

[7] B. I. Abelev, STAR collaboration, Phys. Rev. Lett. 99, 

112301 (2007), and references therein. 
[8] S. S. Adler, et al., PHENIX collaboration, Phys. Rev. 

Lett. 99, 052301 (2007), and references therein. 
[9] B. Alver, et al., PHOBOS collaboration, Phys. Rev. Lett. 

98, 242302 (2007), and references therein. ' 
[10] Bin Zhang, M. Gyulassy, and Che-Ming Ko, Phys. Lett. 

B 455, 45 (1999). 
[11] E. E. Zabrodin, C. Fuchs, L. V. Bravina, and A. Faessler, 

Phys. Lett. B 508, 184 (2001). 
[12] Lie- Wen Chen, V. Greco, Che-Ming Ko, and P. F. Kolb, 

Phys. Lett. B 605, 95 (2005). 
[13] Zhe Xu, C. Greiner, and H. Stoeker, Phys. Rev. Lett. 

101, 082902 (2008). 
[14] Xianglei Zhu, M. Bleicher, and H. Stoeker, Phys. Rev. C 

72, 064911 (2005). 



[15] S. A. Voloshin and Y. Zhang, Z. Phys. C 70, 665 (1996). 
[16] A. M. Poskanzer and S. A. Voloshin, Phys. Rev. C 58, 
1671 (1998). 

[17] Xiao-Mei Li, Bao-Guo Dong, Yu-Liang Yan, Hai-Liang 

Ma, Dai-Mei Zhou and Ben-Hao Sa, to appear in Modern 

Phys. Lett. A; arXiv:nucl-th/0804.3655v3. 
[18] N. Borghini, P. M. Dinh, and J.-Y. Ollitrault, Phys. Rev. 

C 63, 054906 (2001). 
[19] R. S. Bhalerao, N. Borghini, J.-Y. Ollitrault, Nucl. Phys. 

A 727, 373 (2003). 
[20] J. Casalderrey-Solana and U. A. Wiedemann, 

larXiv:09lT 4400vl. 
[21] Ben-Hao Sa, Xiao-Mei Li, Shou-Yang Hu, Shou-Ping Li, 

Jing Feng, and Dai-Mei Zhou, Phys. Rev. C 75, 054912 

(2007); Dai-Mei Zhou, Xiao-Mei Li, Bao-Guo Dong, and 

Ben-Hao Sa, Phys. Lett. B 638, 461 (2006). 
[22] Z. W. Lin, C. M. Ko, B. A. Li, B. Zhang, and S. Pal, 

Phys. Rev. C 72, 064901 (2005). 
[23] T. Sojstrand, S. Mrenna, an d P. Skands, J. Hi gh Energy 

Phys. JHEP05, 026 (2006); |hei>ph/0603175| 
[24] B. L. Combridge, J. Kripfgang, and J. Ranft, Phys. Lett. 

B 70, 234 (1977). 
[25] Ben-Hao Sa and Tai An, Comput. Phys. Commun. 90, 

121 (1995); Tai An and Ben-Hao Sa, Comput. Phys. 

Commun. 116, 353 (1999). 
[26] B. Zhang, Comput. Phys. Commun. 109, 193 (1998). 
[27] B. A. Li and C. M. Ko, Phys. Rev. C 52, 2037 (1995). 
[28] We thank Prof. Lie- Wen Chen for the helps in running 

the AMPT model. 
[29] M. Shimomura, PHENIX collaboration, 

larXiv:0 908.0052v2 . 



